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Abstract 



For the first time the emission of the radiative dissociation continuum of the 
hydrogen molecule (o'^S^ — > b^T,^ electronic transition) is proposed to be used 
as a source of information for the spectroscopic diagnostics of non-equilibrium 
plasmas. The detailed analysis of excitation-deactivation kinetics, rate con- 
stants of various collisional and radiative transitions and fitting procedures 
made it possible to develop two new methods of diagnostics of: (1) the ground 
X^T,'^ state vibrational temperature Tvib from the relative intensity distribu- 
tion, and (2) the rate of electron impact dissociation (d[H2]/c?t)diss from the 
absolute intensity of the continuum. The known method of determination 
of r^ib from relative intensities of Fulcher-a bands was seriously corrected 
and simplified due to the revision of d ^ a transition probabilities and cross 
sections of d <— X electron impact excitation. General considerations are 
illustrated with examples of experiments in pure hydrogen capillary-arc and 
H2+Ar microwave discharges. In pure H2 plasma the values of T^ib obtained 
by two independent methods are in rather good accordance (rvib=3000-5000 
K). In the Il2+Ar microwave plasma it was observed for the first time that 
the shape of the continuum depends on the ratio of the mixture components. 
Absorption measurements of the population of the 3s^3p^4s levels of Ar to- 
gether with certain computer simulations showed that the Ar*— >Il2 excitation 
transfer plays a significant role. In our typical conditions (power flux: 4 W 
cm~^, pressure p = 0.5 mbar, Il2:Ar=l:l) the following values were obtained 
for the microwave discharge: ((i[Il2]/(it)diss ~ 2.5 — 5 • 10^'^ cm~^ s~^. The 
contribution of the excitation transfer is about 10-30% of the total population 
of the state. 
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I. INTRODUCTION 



Experimental and theoretical studies of the emission of the H2 radiative dissociation 
continuum (a'^S^ b^^t electronic transition) have a long and interesting history. They 
were stimulated by molecular spectroscopy and quantum chemistry investigations 
of hydrogen containing plasmas |[^-[T0| and astrophysics ||Tl],|12[. Since the first work of 



Houtermans [13| the — >■ fe^Sj^ transition has been studied as a possible laser system 
(see [p!4| -p!7[] and Refs. therein). It is widely used in UV and VUV light sources [P,p!^-|2T[] . 

In spite of numerous studies, the H2 continuum emission was never used or even proposed 
to be used for a determination of plasma parameters, i.e., for spectroscopic plasma diagnos- 
tics. It seems to be rather strange because, on one hand, diagnostics of hydrogen-containing 
plasmas is not trivial and every opportunity should be tried; on the other hand, emission 
spectroscopy due to its passive nature and instrumental simplicity is of primary interest for 
various studies in plasma physics and industrial applications. Moreover the visible spectrum 
of hydrogen plasma gives us in principle only three qualitatively different sources of infor- 
mation: (1) atomic lines of Balmer series, (2) band (actually multi-line) spectrum of H2, and 
(3) the continuum. First two are successfully used in plasma diagnostics, while the third 
one was never even tried. Previous investigators only measured the continuum intensity and 
sometimes made certain calculations for comparison with an experiment |^,rD|. It should be 
underlined that those calculations are solutions of so-called direct problem of spectroscopy 
(from known plasma parameters to observable spectrum). 

It is obvious that any observable property of plasma is somehow connected with values 
of plasma parameters and therefore in principle may be tried to be used as a diagnostic tool. 
However actually, as a rule, one observable quantity is determined by several different plasma 
parameters. Then the key point is to study the informational content of experimental data 
and to find the proper way of data analysis. The last one is always an inversion problem 
(from an observed spectrum to unknown plasma parameters) having its own specifics (see 
below section |ID . 

The main goal of the present work was to investigate the possibility to use H2 continuum 
intensity for diagnostics of non-equilibrium plasmas. Our decision to make the first steps in 
this direction was based on following a priory considerations. The experimental data on the 
continuum intensity may be considered as two different sources of information: the relative 
spectral distribution and the total absolute intensity of the continuum. 

The first one is directly connected with the population density distribution over vibra- 
tional levels of the a^^g state. In low pressure plasmas the levels are populated predomi- 
nantly by direct electron impact excitation and the deactivation is mainly due to spontaneous 
emission I^JT^ . In this case the relation between the vibrational population density distribu- 
tions in excited and ground electronic states should be obtained within corona-like models. 
So the shape of the continuum should be usable for a determination of the vibrational density 
distribution (or at least the vibrational temperature) of the ground state of H2. 

On the other hand every act of the continuum emission leads to a dissociation of a 
hydrogen molecule. Therefore, an absolute value of the continuum intensity is directly 
connected with the rate of radiative dissociation via the state and also may be used 

for an estimation of the total rate of electron impact dissociation. 



2 



Our intention was to investigate both opportunities and to check them experiment ally.[] 
Therefore the present work includes two closely connected tasks (development of new meth- 
ods and their test in certain real plasma conditions) and two groups of results connected 
with the proposed methods and with non-equilibrium plasma of the discharges under the 
study. 

The first part is based on the analysis and description of: 



1. excitation-deactivation kinetics and derivation of suitable formulas (Sections [11 A| , [11 B| , 

rTD[ , |vg); 



2. rate constants of various collisional and radiative transitions to establish certain sets 
of recommended data (Section C[ ); 

3. experimental technique reproducible by other researchers (Section 

4. various ways of data processing (Section [11 A[ and throughout the paper). 

That made it possible to develop two new methods of diagnostics of: (1) the ground 
X^Sj" state vibrational temperature T^ib from the relative intensity distribution, and (2) 
the rate of electron impact dissociation from the absolute intensity of the continuum. The 
known method of determination of Tyib from relative intensities of Fulcher-a bands was 
seriously corrected and simplified due to the revision of c? — >■ a transition probabilities and 
cross sections of c/ <— X electron impact excitation. 

For tests of the methods we performed the studies of emission spectra of two quite 
different plasma sources: (1) pure hydrogen dc capillary-arc (Section [IV A[) , and (2) H2+Ar 
microwave discharges (Section [IV ij[ ). 

The first one is known to produce bright pure hydrogen spectrum and was already 
investigated spectroscopically P, [I(][ , PT[ , [23[] p[. Therefore it was chosen as an example of a 



simple model system for determination of Tvib by two different methods presented above. 

The second one has been specially designed for basic research on molecular microwave 
discharges mainly based on spectroscopic diagnostic methods p3,|25|. Microwave plasma 



reactors of this planar type have been used for various plasma technological applications 
P3|^H| . Frequently H2 and Ar are used as main components of the feed gas mixture leading 
to an essential interest on H2 dissociation processes and on influences of Ar admixture on 
H2 emission. 

The results are summarized and discussed in the conclusion. 



II. INTENSITIES OF H2 BANDS AND VIBRATIONAL TEMPERATURE 



^Some very preliminary results were already reported in |22]. 



^Our previous attempt to obtain Tvib from Fulcher-a band intensities in 
of insufficient precision of intensity measurements (errors ~ 10%). 



almost failed, because 
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A. Continuum intensity and 
population of upper levels 



The dissociation continuum of the hydrogen molecule is caused by spontaneous transi- 
tions from the upper bound electronic state to the lower repulsive state. The 
spectral distribution /at (A) of the continuum intensity (number of quanta emitted within 
unit range of wavelengths per unit volume and per second in all directions) is related to the 
population densities Nav' of the excited a^S+,t>' vibronic levels, as 

U{X) = Y.Na.'-Af{X), (1) 

v' 

where A - wavelength, v' - vibrational quantum number, 

Nav' ='^Nav'N'-, (2) 



A" 



the total population density of the upper vibronic state, Nav'N' - populations of a S+, t>', A^' 



electronic- vibro- rotational levels (their very small triplet splitting may be neglected), N' - 
rotational quantum number of the total angular momentum excluding electron spin, and 
(A) - spectral distribution of the spontaneous emission transition probability (Einstein 
coefficients for various wavelengths of a^E^,t;' — >■ transitions). We have neglected the 
rather small effect of vibro-rotation interaction and assume the transition probabilities being 
independent on N' . 

For H2 molecule the values of A^^ (A) should be considered as well-known nowadays 
because they have been calculated in several works [ll|,|5|-|7|,p!4| and checked experimentally 



|T|,^, being in rather good agreement 0. We used the results of 0] shown in Table |.0 In 
Fig. |I| the data are visualized for the transitions from a^S+,f'=0-3 levels for A = 200 — 
450 nm. One may see that the probability of transitions from various vibronic levels has 
maxima in different wavelength regions due to oscillatory behavior of the vibrational wave 
functions of the upper state. So the spectral distribution of the continuum intensity /a{,(A) 
may have a certain informational content about the population density distribution over 
vibrational levels of the electronic state. In principle, measurements of the intensity 

/afe(A) may be used for the determination of the populations N^v' in relative or even in 
absolute scale by Eq. ([1|). 

It should be noted that even this very first stage of the data processing belongs to the class 
of so-called reverse (or inversion) problems, which are known to be not well-posed (sometimes 
they are called "ill-posed"). The final result - a certain set of the populations Nav' ~ 
depends (in principle) not only on the values of the actual populations in plasma, but may 
be infiuenced by so-called a priory information: the number of adjusted parameters {Nav'), 



^Those are the only data presented in table form suitable for applications. The data were never 
published in open press, but deposited in VINITI (Institute of Scientific and Technical Information 
of USSR). 
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the number of experimental data points, the experimental errors and even the algorithm of 
data processing. 

In the present work we used for the solution of this and other reverse problems the 
least-square fitting, i.e., the minimization of the functional 

in the multidimensional space of parameters used in a theoretical model. Here 

Y^'^^'^^ai, . . . , a^), Y^^^^' - calculated and experimental data, / - total number of data points, 
m - number of adjusted parameters, cxj - standard deviations of the experimental data. 
The standard deviations of optimal values of the adjusted parameters may be estimated as 
square roots of diagonal terms of the covariance matrix . 

We used a minimization algorithm analogous to a linear regression with the following 
substitutions: Y^^"^^ = /^^^*(Aj) from experiment, = Iab{K) from Eq. (|l]), cr,; = 0.07 • 

J^^''*(Aj) corresponding to 7% random experimental error of the measurement, m = v'^^^^ + l, 
Xi was varied in the range 225 — 400 nm. Experimental errors were estimated by averaging 
of several measurements for the same conditions in plasma. They are mainly caused by 
random noise, non-sufficient stability of the discharges during the experiments and restricted 
reproducibility. The value of v'^^-^ had to be adjusted in such a way that reliable results can 
be obtained from the minimization procedure. The values of A^^ (A) from Table | were 
interpolated by cubic splines. 



B. Excitation-deactivation balance equation 

The populations of Nav' may be sometimes interesting themselves (see the examples 
below). But on the other hand, in non-equilibrium plasmas the rovibronic level populations 
in excited electronic states are connected with those in the ground electronic state by the 
excitation-deactivation balance equation for the electronically excited n', v', N' rovibronic 
levels [0 

^XvN ■ rie ■ ax^N + In'v'N' («,/?, 7, •••) = ^n'^'TV' ( .^d + coH ) ' (4) 
vN W-u'Af' '''n'v'N'/ 

where In'v'N'ic^, • • •) ^ the rate of secondary excitation processes (cascades, recombina- 
tion and so forth); t^?^,^/ - radiative lifetime of the n'v'N' level; t^?^',^/ - effective lifetime 
describing the decay of the n'v'N' level due to quenching in collisions of the excited molecule 
with electrons, atoms and unexcited molecules. The rate coefficient of direct electron impact 
excitation 

= W^)tJ = r at'fie) F{e) da, (5) 

where cf'^^n' i^) ^ corresponding cross section; v,e - velocity and corresponding energy of 
an incident electron; Sth ^ threshold energy of the n', w', N' ^ X,v,N transition; F{e) - the 
electron energy distribution function. 
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For plasma diagnostics this most simple limit case (often called corona or corona-like 
model) is that of low-pressure plasmas with small input power. Under such conditions the 
levels are predominantly excited by direct electron impact excitation and their decay is 
mainly due to spontaneous emission.^ Than the second terms in right and left hand parts 
of Eq. (Q) may be omitted. 

It may be seen from Eqs. (|^) and (|^) that even in the most favorable case for plasma 
diagnostics the populations Nn'v'N' depend not only on the distribution NxvN- Radiative 
transition probabilities, lifetimes and electron impact cross sections for single rovibronic 
levels and transitions should be known as well as the electron energy distribution function. 
Moreover the determination of NxvN from measured N^'v' n' and with the certain system of 
equations for various v' and A^' lead to a reverse problem. 



C. Constants of elementary processes 

Further simplification of the model may be achieved by detailed analysis of radiative 
and collision transition probabilities and lifetimes. We are interested in determination of 
the vibrational temperature Tvib from the continuum intensity (a^S+ b^'^u ) from 
the intensities of Q-branch lines of Fulcher-a bands {d^Hu '^'^^3 )-Q Therefore we have 
to analyze all available data about electron impact excitation and spontaneous decay of the 
upper levels. 

The transition probabilities for d^Il~,v', N'=l — > v", N"=l transitions [Ql lines of 

[v' — v") Fulcher-a bands] were obtained semi-empirically in the framework of the adiabatic 



approximation with corresponding dipole moment obtained in |3o] from experimental data 
about wave numbers, branching ratios in f "-progressions and radiative lifetimes. The results 
for the first seven diagonal {v'=v"=v) bands are presented in Table ^ together with semi- 
empirically predicted [^] and experimental values of radiative lifetimes |^ . 



One may see a noticeable discrepancy between semi-empirical and experimental lifetimes 
of (i^n~, f '=4-6, A^'=l rovibronic levels due to non-adiabatic effects neglected in the adia- 
batic approximation (see also 0,0)- Therefore, only the first four diagonal bands may be 



used for the determination of Tvib. On the other hand this effect did not influence the values 
of the rotational temperatures derived from the populations of d^Il~,v' >3 levels in |Q. It 



may be considered as an indication, that the rate of the additional decay, associated with a 



^An additional population of the a^S^ state due to cascades from e^S^ and d^II.^ states was 



taken into account in [10| and found to be less than 20%. The radiative and collisional quenching 
lifetimes are estimated to be equal at 3 mbar for a^S^ and at 0.75 mbar for d^Il~ states (see 

^The method of Tvib derivation from Pulcher-a bands was originally proposed in It included 
measurements of Q-line intensities, T(e) and gas temperature T. It was based on rather old data 
about radiative lifetimes (see |Q) and the cross sections of the d ■<— X electron impact excitation 
from 1^ . Our analysis show that the data should be revised and the method sufficiently simplified 
(see Eqs. (pi), (10) and the discussion). 
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non-adiabatic coupling of the vibronic levels lying above the H(n=l) + H(n=2) dissociation 
limit, is almost independent on A^' (here n is the principle quantum number). It seems quite 
reasonable, if the perturbation is due to an electron-rotational interaction with a continuum 
(the hypothesis proposed in |^). Because of the lack of information about iV'-dependences 
of the transition probabilities and the lifetimes they will be considered herein as independent 
on the rotational quantum number in accordance with an adiabatic approximation. 

The radiative lifetimes of a^Sj", v' vibronic states of H2 have been studied in numerous 



works both experimentally [p|,^ p!2| , p6| -|38|| and by ah initio calculations All available 

data are collected in Table |TT|. One may see, that: 

(1) Nothing is known about A^'-dependences of the lifetimes, so we again have to assume 
independence on A^'. 



(2) Experimental data are obtained only for low vibronic states. 

(3) The experimental and calculated data are in very good accordance for v' 
for higher v' only ah initio data are available. 



0-2, but 



To be in consistence with the transition probabilities listed in Table | we have used values of 
Tav' marked in Table III as "p.w." - present work. The lifetimes for f ' = — 2 were obtained 
by the integration of the tabulated values of the transition probabilities from Table |. The 
values labeled with an asterisk for f ' = 3 — 6 were obtained by extrapolating the lifetimes 
with the help of a third order polynomial fit of the data from and a constant factor 
obtained from the average value of the ratios with our values for f ' = — 2.^ 

The dependences of ci^II", t>', A^' ^ X^E+,0,A^ electron impact excitation on the incident 



electron energy have been studied in [|32| , ^ 40| . They have a normal form for singlet-triplet 
transitions: a sharp increase from the threshold to the maximum and a rather sharp decrease 
for higher energies. The data of p2| , |39| are in good accordance and in |^T| they were used for 
the extraction of the cross sections for different v' from rate coefficients and F{e) measured in 
0. The relative cross sections o'xoi/^xdi maximum are shown in Table |V] together with 
corresponding ratios of the Franck-Condon factors calculated in |^ . One may see noticeable 



deviation of measured values from those calculated in the Franck-Condon approximation. 



In ||3^ this discrepancy was interpreted as an evidence of a remarkable dependence of the 
scattering amplitude on the internuclear distance (non- Franck-Condon effect). Later the 
data from ^ have been used for the determination of the complete set of the cross sections 
for various f , v' and A^, A^' 



It is important to take into account that the relative cross sections in [^,32,41 



are 



based on radiative transition probabilities and lifetimes obtained in adiabatic approximation 
(almost the same as in Table ^V]) . As we already mentioned above the lifetimes of v' = 
4 — 6 are in contradiction with such extrapolation. New data about t^v'! and from 
Table |I| should be used for the determination of the cross sections from the line intensities 



^This was necessary because for v' >2 one can see in Table | that some quantity of transition 
probability resides below the lower table margin of A = 160 nm. Therefore, the calculation of the 
Tav' from the tabulated (A) values would lead to an overestimation of the lifetimes. 
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measured in M . We made the recalculation and got the data presented in the sixth column 



of Table |^ One may see that they are in rather good agreement with corresponding ratios 
of Franck-Condon factors. It means that the dependence of the scattering amplitude on 
the internuclear distance is negligible in the case of the d^Il~ <— ^^^^ electron impact 



excitation of H2. This is in accordance with recent ab initio calculations ||43[| . 



Nothing is known, so far as we know, about the cross sections for the a^S+,f',A^' ^ 
X^J^t ,v, N excitation. Therefore, for both d^Il~ <— X^Sj" and a^Sj" ^ X^J^f transitions 
we assume that the rate coefficients of the electron impact excitation have non-zero values 
only for N' = N rovibronic transitions. They are considered as independent on N and 
proportional to corresponding Franck-Condon factors, i.e., 

"x^/ ~ « Qx'', when AN' = N'-N = (6) 
and ai^Jf = for AN' ^ (7) 



thus neglecting the momentum transfer in electron impact excitation [32,44| and the rather 



small dependence of the rate coefficients on F{e). We used the only values of Franck-Condon 
factors for d^Il~ ^ . 
(Tables [V| and |Vip.[] 



factors for d^U,, X^H^ and a^S^ <— X^H'^ transitions which had been calculated in |^ 



D. Determination of vibrational temperature 

Vibrational and rotational energies are not well separated in the ground electronic state 
of H2 because of the small mass of the nuclei. Just to have a certain ability to characterize 
the distribution over vibrational levels we introduce the vibrational temperature by 



^ V/V v- P2]-exp(-^).(2iV + l).exp(-A^) exp ^ 



N N 



Eexp(-fg-) E(2iV+l)-exp(-^) ^^vib) ' 

(8) 

where [H2] - total concentration of molecules, AExv = Exvq — Exoo, ^ExvN = ExvN — Exvo 
- vibrational and rotational energy differences, T^ib and T - vibrational and translational 
temperatures, 'P(Tvib) - the vibrational partition function. So we assume that in all vibra- 
tional levels of the X^T^^ state the populations of the rotational levels are in Boltzmann 
equilibrium with the gas temperature T. 

Taking into account all the assumptions discussed above, the balance equation (^) may 
be written for the populations of (i^n~,f' vibronic states after summation over rotational 
levels in the following form 



^Although those data should in principle be available by order from Russia, we decided to include 
them into the article for completeness of the data set on elementary processes involved and for 
convenience of potential users of our methods. 
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One may see that the determination of T^n, may be achieved by numerical solution of 
the system of equations for several levels v' with experimental data about N^v' derived 
from the intensities of the Fulcher-a Q-branch lines of diagonal (f ' = v") bands: 

Tdv'N 

^^^' = E#^- (10) 

This method of the Tvib determination is based on sufficiently simple kinetics model and 
certain set of constants presented in Tables || and |VI|. It is even more simple then the 
predecessor because it needs only intensity measurements. 

In the framework of our model the continuum intensity can be written as 

1 '^max / /\ fP \ 

UX,T.i^) « :^7;^E^r'(A) ■ W • E Qt ■ exp [-^). (H) 

To be able to work with relative intensities we may introduce the normalized relative 
intensity 

<^b(A, Ao, Tvib) = /ab(A, rvib)//a6(Ao, Tvib), (12) 

which is equal to unity for A = Aq. 

The normalized intensities (|12D may be tabulated for various values of Tyjb and then 
may be used for a determination of the vibrational temperature by least squares fitting 
(0) of experimental and calculated continuum intensity distributions. The results of our 
calculations for Aq = 355.5 nm and Tvib = 0, 3000 and 5000 K are shown in Fig. ^ One 
may see that a variation of vibrational temperature leads to certain change in the shape 
of the continuum intensity distribution. Most sensitive are the wavelength regions near 
A ~ 200 and 300 nm. The first one is out of our range of observation and cannot be used 
for diagnostics herein. Nevertheless the fact should be considered as one of the results of 
our present work. It may be recommended for VUV spectroscopy if the problems of the 
sensitivity calibration in this difficult region would be solved somehow (with the use of 
synchrotron radiation, for example). The changes near A ~ 300 nm are less pronounced, 
therefore they may be used in the following way. 

When the vibrational temperature is sufficiently low and 

Qx>j:Qx:-e-^, (13) 



It's also an inversion problem, but its solution is rather stable because it is quite easy to measure 
intensities of the first four bands (u' = — 3) . That means at least three relative populations ( p!o|) , 
and therefore three equations (^) can be used for derivation of one unknown (T^ib). 
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then the excitation from vibronic levels with v > may be neglected and 



E (A) ■ r,,, ■ 

Jab{\, Ao, 0) = ^ (14) 

J2 (Ao) ■ Tav' ■ QxO 

v'=0 

represents the shape of the continuum for T^ib — > 0. It does not depend on discharge 



conditions and may be easily calculated because molecular constants in Eq. (14) are known. 
The results of such calculation for Ao = 355.5 nm and t>^^^ = 6 are presented in the 
last column of Table |. One may see that in the wavelength range of observation, A = 
225 — 400 nm, the continuum intensity shows monotonic increase towards short wavelengths. 
The ratio 



J'abi-^^ Ao, Tvib) — Jabi^, Ao, T^ih)/ Jabi^, Aq, 0) (15) 

may also be calculated for various values of the vibrational temperature and then be used as a 
nomogram for the determination of T^ib from measured relative continuum intensity without 
solving the inverse problem by a minimization of Eq. (^. The results of such calculations 
for a limited number of T^ib = — 6000 K are shown in Fig. 0. One may see that the ratio 
(p!5|) has some sensitivity to the vibrational temperature especially near 300 nm. For this 
particular wavelength we may organize a graph of ^^^(A = 300 nm, Aq = 355.5 nm, T^ib) 
versus T^ib presented in Fig. ^ where the sensitivity is even more clear. This dependence 
may be used as a nomogram too, being especially useful for in situ control of plasmachemical 
processes. 



III. EXPERIMENTAL SETUP 



As it was already mentioned in the introduction the experiments in pure hydrogen and 
H2+Ar plasmas were carried out in two different gas discharges. 

The first one is a hot cathode arc with 2 mm diameter constriction analogous to those 
used in ||19| , |20[| . The discharge device was filled with 8 mbar of spectrally pure hydrogen. 



The range of the discharge current / = 10 — 500 mA allowed us to achieve rather high current 
densities j ~ 0.3 — 16 A cm~^. 

The second plasma source was a H2+Ar microwave discharge excited by a planar mi- 



crowave applicator similar to that described in 1 24 1 . This plasma reactor has been specially 



designed for basic research on molecular microwave discharges mainly utilizing spectroscopic 
diagnostic methods |^ . The discharge configuration has the advantage of being well suited 
for end-on spectroscopic observations, because considerable homogeneity can be achieved to 
a certain extent, while side-on observations are also possible. Microwave plasma reactors of 
this planar type have been used for diamond deposition p6| , surface corrosion protection by 
deposition of organo- silicon compounds and surface cleaning |^S[. The inner dimensions 
of the discharge vessel are 15x21x120 cm^ and the area of the quartz microwave windows 
was 8x40 cm^. The microwave power was fed in by a generator (SAIREM, GMP 12 KE/D, 
2.45 GHz). We were able to measure the total input (W = 1.5 kW) and reflected power. 
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the latter was in our conditions always negligible (< 4%). The input power flux in the plane 
of the applicator windows was roughly estimated to about 4 W cm~^ assuming 20% heat 
loss to coolant and environment. The plasma length in end-on direction was 15 ± 3 cm, 
mainly depending on the discharge pressure, and 8 ± 1 cm in side-on observation, while 
for absorption spectroscopy the end-on direction was preferred. Gas flow controllers and 
a butterfly valve in the gas exhaust allow to have an independent control of gas inflow (y9 
= 100 seem) and pressure (p = 0.5 — 5 mbar). The purity of gases was 99.996%. H2 and 
Ar are used as components of the feed gas mixture, since they are not only of interest for 
basic research but of importance for various applications as well. Our main interests were 
focused on studies of: (1) influence of Ar on H2 emission, and (2) dissociation of molecular 
hydrogen. 

The central part of the plasma along the axis of the capillary discharge was focused 
by a quartz achromatic lens onto the entrance slit of a monochromator. In the microwave 
plasma the line of sight was 12 mm below the microwave windows of the applicator. For 
measurements of the continuum intensity we used a 0.5 m Czerny- Turner monochromator 
(Acton Research Corporation Spectra Pro-500) with gratings of 600 and 1200 grooves per 
mm in the flrst order. The line intensities of the Fulcher-a bands however were measured 
using a high resolution 1 m double monochromator of Czerny- Turner type with 1800 grooves 
per mm (Jobin Yvon RAMANOR UIOOO) in the flrst order.[] The output light was detected 
by a CCD matrix detector (TE-cooled, 512x512 pixel, size of image zone: 9.7x9.7 mm^) 
of the Optical Multichannel Analyzer PARC OMA IV connected with a computer. Special 
software made it possible to collect and analyze the recorded spectra. 

The wavelength calibration was done by means of known Hg, Ar and II2 spectral lines. 
The relative spectral sensitivity of the spectroscopic system in the range A = 225 — 450 nm 
has been obtained using a deuterium lamp (Heraeus Noble Light D0651MJ) calibrated 
by the Physikalisch-Technische Bundesanstalt (Braunschweig). The absolute value of the 
sensitivity for A = 350 nm was determined with a tungsten ribbon lamp. Its current- 
temperature calibration has been performed by the Russian State Institute of Standards 
(St.-Petersburg). We used two different currents corresponding to the absolute temperatures 
of tungsten T = 2637 and 2763 K. This gave us usually a spread in the values of the sensitivity 
(at A = 350 nm) of less than 1%. The emissivity coefficient of the tungsten ribbon for these 
temperatures was taken from [^] and used for the calculation of the spectral distribution 
of emission intensity. 



^ The role of atomic hydrogen is usually considered as very important for plasma-surface interac- 
tion processes [^,46], but it was not investigated in the planar microwave discharge. 



^*^We started our measurements with the grating of 600 grooves per mm (blazed for 300 nm) 
because of higher sensitivity [see Fig. |5|(a)] important in the case of microwave discharge with its 
low intensities [pay attention to ordinate scales of Fig. |5|(c) and P(a,b)]. Later we had to switch 
to the grating with 1200 grooves per mm (blazed for 250 nm) to get better relative sensitivity for 
shorter wavelengths [see Fig. |5|(a)] and to select the continuum intensity from that of the band 
spectra of H2 in more accurate way. 
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Some typical experimental records of the spectra of a standard deuterium lamp (a), pure 
hydrogen capillary-arc discharge (b) and H2+Ar plasma of the microwave discharge (c) are 
presented in Fig. ^ The non-monotonic structure in the range A = 250 — 350 nm is connected 
with the spectral distribution of the sensitivity of our spectrometer. In the long-wavelength 
end of the first two records the band (actually multi-line) spectra of D2 and H2 molecules 
may be seen. One may see also that in spite of the gas flow mode our microwave plasma 
is not spectrally pure. It contains some impurities: oxygen (283 nm), OH-radical (bands at 
283, 310 — 320 nm), NH-radical (336 nm) and N2 molecule (bands of the 2nd positive system 
280 — 400 nm mixed with H2 bands). It should be noted that some lines are overexposed in 
the record shown in Fig. |^(c). The data of the continuum intensity, except those in Figs. 
^ and ^, were manually extracted from the measured intensity distribution. 

As can be seen from previous description, the spectroscopic part of our experimental 
setup is certain combination of components commercially available nowadays. Therefore 
it may be easily reproduced. Experimental errors are mainly caused not by the detection 
system but by the random noise of plasma sources used. Special attention should be paid 
to selection of the continuum intensity from that of multi-line spectrum of H2. 

The results of our experiments will be presented together with corresponding analysis in 
the following sections. Perhaps most interesting of them are the first observations of: (1) 
small changes in the shape of the continuum in pure hydrogen plasma caused by vibrational 
excitation in the ground electronic state of II2 molecule,|^ (2) influence of Ar concentration 
on the shape of the continuum in H2+Ar plasma of microwave discharge. 



IV. RESULTS AND DISCUSSION 
A. Pure hydrogen dc capillary-arc discharge 

Relative intensities of the continuum and of the first Q-branch lines Ql to Q5 of Fulcher- 
a bands (f '=f"=0-5) were measured in the discharge current range mentioned above. The 
intensities were found to be proportional to discharge current if one takes into account a 
decrease of the density of molecules due to warming of the gas inside closed lamp by the 
hot cathode and discharge current. As an example the signals of CCD detector proportional 
to the continuum intensity are shown in Fig. |^ as a function of the discharge current. The 
linearity observed may be considered as an argument in favor of direct electron impact 
excitation of the levels under the study. 

The shape of the continuum intensity distribution was observed to be independent on the 
discharge current for currents i < 100 mA within experimental errors. A typical result used 
in further analysis is presented in Fig. |^. As one may see, measured continuum intensity 
is close to that calculated for T^ib 0. The low wavelength margin is limited by a sharp 



^^That was achieved due to sufficient improvement of experimental technique in comparison with 
previous investigations ||8-1C| where the changes have been masked by random noise of discharge 
and large systematic errors. 



12 



decrease of the sensitivity of our spectrometric system. For wavelengths longer than 360 nm 
the overlap of the continuum with band spectra of H2 becomes more and more pronounced. 

The first question to answer is how many populations of upper vibronic states can be 
obtained from the intensity distribution measured in the wavelength range of observation. To 
find the optimum number of adjusted parameters Nav' we made the x^-minimization [Eq. (^] 
with various m = (f ^ax + The results of the calculations, made for the experimental data 
of Fig. ^ are presented in Table |Vli| . shows a monotonic decrease and becomes lower 



than unity for f^^^ > 3. Another characteristic, presented in the table, is 



Pmax = max(pj) with pi = ^l-l/^i, (16) 
where ^^ = ZuZr^' and = (17) 

The magnitudes pi are characteristics of the precision of the solution. Its zero value corre- 
sponds to the limit case of total independence of the parameters, while Pi = 1 means that the 
parameter may be expressed as a linear combination of the other parameters and therefore 



the number of adjusted parameters should be decreased. It may be seen from Table |VII 
that in our case pmax is growing with the increase of m; pmax is close to unity for v'^^^ >3. 
The maximum value of the relative standard deviation of the parameters {ANav' /Nav')ma,x 
(in our case it is usually observed for v'=v'^^^) is small and almost independent on v'^^^ for 
■^max = — 3, then it jumps to about 0.83 at f ' = 4 and becomes much bigger for v' = 5, 6. 
It means that the standard deviation of Na4 derived from our spectrum is almost the same 
as its optimal value. 

Thus in our case a determination of the populations of more than four vibronic levels is 
meaningless. The reason is very simple and may be seen from Table |V111| where the relative 



populations Nav'/Nao (obtained for v'^^^^ = 3) are shown together with relative intensities of 
various v' b^^u transitions in the spectral range under the study (A = 225—400 nm) 



UX) dX 



= Xr .. :r ■ (18) 



From the data of Table |V111| it becomes clear that only 4% of the detected emission quanta 



originate from levels with v' > 3. The populations Nav' observed in the pure hydrogen 
capillary-arc discharge for currents i = 10 — 100 mA are quite close to those calculated in 
the framework of our simple model with T^ib — > 0. It means that the populations of the 
X^S^, V vibronic levels with f > are so small that they do not affect relative populations 
of a^S^,f' levels and therefore relative continuum intensity distribution. 

For currents higher than 50 — 100 mA we observed some systematic changes in the 
shape of the continuum intensity distributions. As an example our experimental data on 
the normalized intensities J^j,(300, 355.5, T^ib) from Eq. (p!5D are presented for three values 
of the discharge current in Fig. ^ together with the curves calculated with our model for 
various values of T^ih- One may see that two main predictions of our model - the maxima 
near 300 nm and crossing with J^^ = 1 near 250 nm - are in qualitatively good accordance 
with the experiment. The experimental curves are a bit shifted to shorter wavelengths. 
To treat the data quantitatively one should take into account that the differences between 
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experimental and theoretical curves are comparable with errors of both the experiment and 
modeling: 



(1) The actual random noise of our measurements is about two times bigger than that 
shown in Fig. ^ due to smoothing during calibration procedure. It was made for 
better visibility of the trend. 

(2) The normalized intensity J^^ is the ratio of the measured intensities, so its errors are 
twice more than that of measured intensities. 

(3) The systematic errors of the measurement may be even bigger because errors of the 
sensitivity calibration should be included, while the errors of the relative energy cali- 
bration of the deuterium lamp are about 4%. 

(4) Errors of the calculations are not so easy to estimate, but most likely they may be 
of the same order of magnitude. They are caused by the approximations made in 
the formulation of our very simple excitation-deactivation model [Eq. (Q)] and by the 
uncertainties of the rate coefficients (see Section p.1 Cj ). 



Nevertheless, we used the values of maximum deviation of experimental curves from J^^ = 1 
level for estimation of T^ib- The results are shown in Fig. ^ together with values of Tyib 
obtained from Fulcher-a band intensities by a system of Eqs. (|^) for v' = — 3.0 One may 
see that both spectroscopic methods of Tyib determination developed in the present work 
are in good accordance. 

The ground state rotational temperature was obtained from the intensity distributions 
in Fulcher-a Q branches by the method proposed in PSI. In our conditions the rotational 



and translational (gas) temperature are known to be in rather good accordance p2| , |30| , |34 
The results are also shown in Fig. |^ together with earlier data obtained in |^ . One may see 
that our measurements show remarkable difference between Tvib and T analogous to that 
observed in other conditions by CARS |^ |. 



For currents higher than 300 mA the continuum intensity distribution became again 
insensitive to the current variation. This is most probably caused by large values of gas 
temperature and concentration of atomic hydrogen that is favorable for the acceleration of 
the VT-relaxation. On the other hand this behavior can also be explained by our model, 
since there is a certain upper limit above which the continuum intensity distribution is no 
longer sensitive to an increase of vibrational temperature (see Fig. 

B. Microwave discharge in H2-|-Ar mixture 

The results of our measurements show that the spectral intensity distribution of the 
continuum in the H2+Ar microwave discharge has a noticeably different shape compared 
to that in pure hydrogen. Moreover, the shape was found to depend on the ratio of the 



^^Intensities of first five lines Q1-Q5 were used in each band. 
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components in the mixture ([Ar]:[H2]), which was varied from (19:1) to (1:4) under the same 
total pressure and power input. We observed a relative increase of the continuum emission 
in the wavelength region near 300 nm. As an illustration some typical results for (4:1) and 
(1:4) mixtures are shown in Fig. One may see that the additional emission is located in 
the same wavelength interval as a^Ti^ ,v'=0—>- transition of H2 (see Fig. |1]) and ArH* 
(see below). Note that measurements in between exhibit a monotonic dependence of the 
increased emission in the wavelength range around 300 nm. An analysis of various possible 
explanations of the observed phenomenon lead us to the assumption that in our conditions 
an excitation transfer from excited Ar atoms to the hydrogen molecules might take place. 
This experimental observation leads to two important consequences: 

(1) The shape of the continuum in H2+Ar plasma is affected by two different excitation 
processes, therefore our model based on one of them (direct electron impact excitation) 
can not be used for the determination of Tvib.Q 

(2) To check the hypothesis of excitation transfer in our plasma we have to accompany 
our emission intensities with absorption measurements of populations of long-living 
states of Argon atoms (see Section |1V B 2| ). 



1. Ar*^H2 excitation transfer, prehistory 

Excitation transfer in collisions of metastable (^Po2) resonant (^P^' and ^P^') Ar 
atoms with hydrogen molecules is known since the pioneer work of Lyman and was 
investigated so far in plasma and beam experiments mainly as: (1) quenching of excited Ar 
and (2) additional excitation of VUV Lyman and Werner bands (-B^S+, C^Tlu "^^^g ) 
of H2. The first attempts to study the role of this effect in the emission of the hydrogen 
dissociation continuum in H2+Ar plasmas appeared quite recently [ P^ , ^T| -|55 



The excitation of a S+ state of H2 due to the excitation transfer from excited Ar* atoms 
to the ground state H2 molecule may go in two different ways: 

Ar*(4si'3p}) + H2(XiS+ t;,Ar) + To ^ (ArH2)* ^ (19) 
H2*(a3S+ v', N') + Ar(iSo) ^ 2 H(l2Si/2) + Ar(iSo) + Ti + /^//(Hs) (|19la) 

ArH*(A2S, v', N') + H(l2Si/2) ^ 2 H(l2Si/2) + Ar(^So) + T2 + /iz/(ArH). (||6) 

Here (ArH2)* denotes a temporary excited state of the system during the interaction, Tq and 
Ti, T2 are the initial and final kinetic energies of reacting particles in the mass center frame. 
The branching ratio between these two output channels of the reaction should depend on 
the initial excited state of the Ar atom and the vibro-rotational state of the molecule as well 



^The further analysis shows that this loss is well compensated by the opportunity to get informa- 
tion about the mechanism of excitation of the continuum emission and to account for the excitation 
transfer in determination of the total dissociation rate. 
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as on the collision energy Tq. As far as we know it is not established up to now in spite of 
a certain attempt in crossed beam experiments [Q. 

It may be shown, that if one wants to consider the excitation transfer possibility for ^'^P} 
argon atoms with various J values from the energy point of view, the existence of vibro- 
rotational and translational energies of H2 molecules in plasma should be taken into account 
|53| . Then in our conditions the observed effects are most likely due to the collisions of the 
resonant ^P^ and metastable ^Pg argon atoms with the X^E^,f=0 hydrogen molecules. 



2. 3s^3p^4^ level populations of Ar 

To clear up the situation we measured the population densities of the Ar 3s^3p^4s (^Pg 1 2 
and ^P°) levels by an ordinary self-absorption method with one mirror behind the discharge 
vessel g. 

Neglecting Stark broadening of the spectral lines the population of the initial state of 
k <— i absorption transition may be presented as 



where 

/o = ^ = 2.64 ■ 10-2 cm2 s'S (21) 



vre^ 



nipC 



fik - the oscillator strength for the k ^ i transition in units of /o, Auo - the Doppler 
half-width depending on the gas temperature. 

The absorption coefficient kq was obtained by numerical solutions of the following non- 
linear equation [^ : 



\ _ S{2koI) 



:^ + ^)-^rn/l0^ (22) 



S{koI) 

where Im and Jq - the line intensities measured with and without the mirror, / - length of 
plasma column along the axis of observation, r - the reflection coefficient of the mirror and 
S{x) - the Ladenburgh-Levy function calculated by formulas from 

The effective reflection coefficient has been determined experimentally by the measure- 
ments of several Ar lines (603.1, 703.0, 714.7, 789.1, 860.6, 862.0, and 876.1 nm) which are 
free of self-absorption. It was found to be r = 0.52 ± 0.02 and independent on A. We used 
the values of the oscillator strengths from ||56[| . 

The intensities of two or three different spectral lines were used for the population den- 
sity determination of every 3s23p^4s level, namely: 

3p° - 763.511, 801.479, 811.531 nm 
3p° - 751.465, 738.398, 810.369 nm 
3pg - 794.815, 866.794 nm 
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^Fl - 826.452, 840.821 nm. 



The conditions were the same as those used in our measurements of the continuum intensity. 

Typical results are shown in Fig. ^ The scaled populations are absolute values of the level 
population densities divided by their statistical weights (2 J + 1). They have the meaning 
of the average populations of the multiplet structure sublevels. In equilibrium conditions 
they should be almost equal to each other because the levels have almost the same energy - 
11.55, 11.62, 11.72, and 11.83 eV for 3p°, 3p°, 3po^ and ^P^. One may see from Fig. | that 
in our conditions: 

1. The populations of various sublevels of all resonant and metastable states coincide 
within the experimental errors. It is a direct manifestation of long effective lifetimes 
of the resonant levels (due to self-absorption of the resonant radiation) and a high 
enough rate of collisional transitions between the 3s^3p^4s levels. So the metastable 
and resonant levels of Ar are actually mixed in the plasma under the study. 

2. The sub level populations show the monotonic decrease with total pressure which may 
be connected with: a decrease of microwave power applied to plasma in the volume 
of observation, a decrease in the rate coefficient of electron impact excitation (due to 
certain changes of electron velocity distribution) and an increase of collisional quench- 
ing. 

3. The total population of all 12 sublevels may reach 10^^ cm~"^ at low pressures. The 
rate coefficients of the Ar*— excitation transfer are approximately 10~^° cm^ s~^ 
[0. Thus the rate of the continuum excitation in the reaction is about 10^^ cm~^ s~^. 



This is comparable with the observed total continuum emission (see Table fX\} and 
should be taken into account together with direct electron impact excitation of the 
a^S^ state of H2. 



3. Shape of the continuum in H2 +Ar plasma 

Let us assume the spectral distribution of the continuum intensity as the sum of three 
independent terms corresponding to three different excitation mechanisms: 

= A[cJab{\ 0) + C2lt\\) + C3/ArH(A)], (23) 

where A - absolute scale normalization constant equal to total continuum intensity in the 
range of observation A = 225 — 400 nm, /ab(A, 0) - spectral distribution of the continuum 
radiation caused by the electron impact excitation of H2 from lowest X^S+,f=0 vibronic 
ground state [Eq. ([T^ ], /^°(A) - the same due to the reaction (|T9|a) for t>' = only, /ArH(A) 
- the spectral distribution corresponding to the reaction (|19]b) calculated in i, a - the 
coefficients representing relative contributions of the components iJZiCi = 1). The spectral 
distributions /ab(A,0), and /ArH(A) where normalized for unit area in the range of 

wavelengths A = 225 — 400 nm. 

Then the spectral distributions of the continuum measured in the microwave plasma 
under various conditions were fitted by the function ( ^3]) according to Eq. (^. Typical 



17 



results are presented in Table |^for five selected mixtures of Ar+H2 from (1:4) up to (19:1). 
The first set of x^-minimizations marked as (ci, 02,03 7^ 0) was made without any additional 
preconditions. One may see that the observed shape of the continuum may be described by 
three adjusted parameters of the expansion (p3D with rather good precision (x^ = 1 — 2). 
Moreover the fittings show the negligible contribution of the reaction (p!9|b), which cannot be 
considered as significant because the values of C3 are comparable with standard deviations 
of their determination. Figure 110 illustrates typical results of such fitting. The contribution 
of the reaction (|19|b) is too small to be visible in the figure. 

The second set of fittings has been performed under the additional condition C3 = 0. This 
version of the model provides almost the same precision in the description of the observed 
spectra taking into account only the first two mechanisms included into Eq. (pSj) . 

In contrast the third set of fittings has been carried out for C2 = thus neglecting a 
contribution of the reaction (|l^a). Again we got almost the same quality of fitting. These 
results are the direct consequence of the similarity of spectral distributions of the continuum 
emission caused by the reactions ([T9|a) for v' = and (|19]b). 

It should be noted, that ab initio potential curves of ArH* where never checked experi- 
mentally and the spectral distribution of the ArH continuum has been calculated in 0] in the 
rough approximation of harmonic oscillator. The precision of such calculations cannot be 
high nowadays. Taking this into account we may assume that the insignificant contribution 
of the reaction (|l9|b) obtained in the first set of fittings is an accidental numerical result. 
Most probably it is without any physical meaning because two other sets [based on directly 
opposite assumptions about the branching ratio between reactions (|T9|a) and (|T9|b)] gave us 
almost the same values. 

Therefore, we have to conclude that the precision and the wavelength range of our 
measurements together with the lack of dependable ArH* transition probabilities do not 
allow us to distinguish contributions of the reactions ([T9|a) and (|l9|b). On the other hand, 
there is no doubt that in our conditions the role of the Ar*— >-H2 excitation transfer is well 
noticeable (10-30%) both in excitation and dissociation rates (see below). 



C. Estimation of the rate of electron impact dissociation 

As it was already mentioned in the introduction, the absolute value of the intensity of 
the continuum is directly connected with the rate of the radiative dissociation via sponta- 
neous emission to the repulsive state. If self-absorption in the plasma volume may be 
neglected, then the rate of radiative dissociation due to transitions is equal 

to the total continuum intensity (expressed in number of photons per cm^ in a second) 
integrated over the entire wavelength range 




(24) 

The wavelength range where we were able to measure the intensity was limited (see Fig. 0, 
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H, H). So we detected only a part of the continuum emission.0 The rest should be obtained 
by some extrapolation. In pure hydrogen discharge plasma the vibrational temperature 
is usually rather low and the Eqs. (|1|) and (|12D may be used for the extrapolation. The 
calculations show that in our case: (1) The dependence of the total intensity (area under 
the curve) on the v'^^^ is insignificant after v'^^^ became high enough. (2) The contribution 
of the undetected part of the spectrum is rather significant (about 60% in the example shown 
in Fig. I). 

The other, more general way of treating of spectroscopic data is a numerical solution of 
the reverse problem by x^-minimization (^). The relative or absolute populations Nav' thus 
obtained may be used for the determination of the relative or absolute radiative dissociation 
rate. Radiative lifetimes are known, so an integration over wavelengths is not necessary and 
the procedure is reduced to only few summations in Eq. (p^). 

In low-pressure gas discharges the upper levels of the continuum transition are predom- 
inantly populated by *— X^E+ electron impact excitation [^|Tni- A dissociation of 
hydrogen molecules is also mainly caused by electron impact excitation of various inter- 
mediate excited electronic states, bound and repulsive ones. The relative contribution of 
various channels to the total electron impact dissociation cross section cTdiss(£^) was analyzed 
WJ\ , later supplemented with an application of the close-coupling calculations of the cross- 



m 



sections |E8[. It was found that singlet-triplet transitions play the major role, <— X^S 



g 



and a^E+ ^ -^^^^ being most prominent. The experimental data for crdiss(^) are available 



from p9 |. 



All these cross-sections seem to have a rather similar shape as a function of the collision 
energy e, except in the small region near the threshold, because an excitation of various 



electronic states obviously has different start-up energies. In |22| we proposed to use average 
values of ratios between the cross sections (determined far enough from threshold region) 
for the estimation of the rate of the radiationless ^ -^^^t Process ((i[H2]/(it)f,^js: 

and of the total dissociation rate ((i[H2]/(it)diss by the multiplication of the measured total 
continuum intensity with certain coefficients (4 and 15 correspondingly). Actually, the 
coefficients depend on the shape of the electron energy distribution function and the total 
electron impact dissociation rate 




diss 




(25) 



Generally, the rate coefficients in ( p5|) should be calculated with actual F{e) by Eq. (||) P,p!0 



If the distribution function is not available like in our present study, only rough estimations 
are possible. It is quite obvious that, when e-dependences of the cross-sections have similar 
shapes, the ratio of the rate coefficients is less sensitive to the near-threshold region for F{e) 
with a high content of fast enough electrons. To check the tendency we calculated the ratios 
for a Maxwellian distribution with an electron temperature of Te = 1 — 10 eV and the cross 



Measurements in VUV are possible in principle, but they are not typical in most gas discharge 
and plasma technology applications due to their technical complexity and the difficulties in the 
sensitivity calibration and the spatial resolution. 
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sections from |l58|j59| . The results are presented in Fig. |TT[ One may see that both ratios show 
sharp decrease for < 5 eV becoming almost independent from for Te > 5 eV. It is quite 
typical for a low-pressure hydrogen discharge plasma to have a strongly non-Maxwellian F{e) 
with a high density of fast electrons having energies sufficiently higher than the thresholds 
(see i.e. [pO|). So we may use the asymptotes of Fig. O for the estimations. Numerical 



values of the coefficients are 

^^^4.9 and ^^^16. (26) 

In the pure hydrogen capillary-arc discharge we were not able to carry out absolute 
measurements because the plasma is inhomogeneous along the axis of observation and the 
effective length of the emitting plasma column is uncertain |l60 |. 

However, in the H2-|-Ar microwave plasma the length of emitting column (8 ± 1) cm is 
well defined and absolute continuum intensities were measured in a certain range of discharge 
conditions. But the Eqs. ( P^ and (|2^) are not valid any more, because Ar*— excitation 
transfer plays an important role. The continuum intensity may be expressed as (p3|) and the 
total dissociation rate is a sum of the terms corresponding to the electron impact and the 
excitation transfer. It can be shown that in this case the total dissociation rate 

fd[}l2]\ 



dt 




(27) 

total \ / diss \ ""^ / Ar* 

oo /"OO _ /■oo 

(28) 
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Spectral distributions of the intensities in the second and third terms are mainly located in 
the wavelength range of our observations (about 90%). Only the first one, connected with 
a^E+,f' ^ X^S+,t>=0 excitation needs the extrapolation described above. Therefore, the 
total dissociation rate in our conditions may be calculated from the measured continuum 
intensity. The results are shown in Fig. |T2](a,b) for various pressures and ratios of the 
components in the mixture. They may be interpreted only qualitatively. 

Figure 0(a) shows that the dissociation rate in the volume of observation monotonically 
decreases with increasing pressure. The effect is analogous to that of the excited argon 
populations (see Fig. ^ and corresponding discussion). It is caused by the mechanisms 
already mentioned above, mainly by the decrease of microwave power coming to the volume 
of observation. 

The observed changes in the dissociation rate caused by the variation of mixture compo- 
nents under constant pressure [Fig. |l5(b)] are not trivial. In the absence of secondary effects 
one should expect linear dependence, whereas our measurements show a non-linear decrease 
of {d[B.2]/dt)diss towards smaller values of [H2]. Most probably the effect is connected with 
very simple reasons, which can be of importance for industrial systems. When the content of 
hydrogen in the input gas mixture is high enough, the concentration of hydrogen in plasma 
is proportional to that in the input gas mixture. If the hydrogen content in the feed gas is 
reduced, then the hydrogen density in plasma is mainly influenced by resorption of hydrogen 
and water from the rather large metal surfaces.^ If it would be necessary to achieve small 



The concentration of desorbed water vapor has been measured by IR absorption spectroscopy 
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concentrations of hydrogen in plasma, then heating and pumping of the reactor must be 
done for a very long time. 

On the other hand the data presented in Fig. p!2|(b) show that in our conditions the 
Ar*— i>H2 excitation transfer plays noticeable role even in the total rate of hydrogen disso- 
ciation. This effect should be taken into account in kinetics modeling of atomic hydrogen 
concentration in plasmachemical systems where Ar is used for transportation of reagents as 
a buffer. Our observation shows that it is not just a buffer. 

V. CONCLUSION 

In this paper for the first time the H2 radiative dissociation continuum was used as a 
source of information about parameters of non-equilibrium plasma. Two new methods of 
spectroscopic diagnostics of hydrogen containing plasma have been developed: (1) deter- 
mination of the vibrational temperature of the H2 ground state from the relative intensity 
distribution and (2) derivation of the rate of electron impact dissociation from the absolute 
intensity of the continuum. The development of these methods has been based on a detailed 
analysis of the excitation-deactivation kinetics, rate constants of various collisional and ra- 
diative transitions as well as the way of data processing. The known method of vibrational 
temperature determination [^] using H2 emission line intensities of Fulcher-o; bands was sig- 
nificantly improved and simplified. The potential of the new methods for plasma diagnostics 
was demonstrated at pure H2 capillary-arc and H2-|-Ar microwave discharges. 

In pure hydrogen discharge plasma it was observed for the first time that for high enough 
current densities j > 3 A cm~^ the shape of the continuum intensity distribution is influenced 
by the vibrational temperature of the ground electronic state. Two independent 

spectroscopic methods gave almost the same values of Tvib ~ 2000 — 5000 K. The difference 
between vibrational and gas temperatures is in accordance with previous observations ^ 



In the H2+Ar microwave plasma for the first time the shape of the H2 dissociation 
continuum was found to depend on the mixture components, significantly influenced by 
Ar* — s>H2 excitation transfer processes. Absorption measurements of the population of the 
3s^3p^4s levels of Ar together with certain computer simulation has been successfully used 
for verification of the proposed mechanism. Also for the first time the rates of radiative 
and total electron impact dissociation were obtained for various gas pressures and relative 
concentrations of the mixture components. 

Our attempt to study the informational content of the intensity of the hydrogen disso- 
ciation continuum and to check its applicability for plasma diagnostics has been a priory 
limited in two different respects. 

1) We intended to have an experimental technique as simple as possible and, in par- 
ticular, to stay on the basis of emission spectroscopy only. Then the methods would be 
especially useful in various industrial applications including rather easy in situ control of 
plasma processing. The approach has in principle a backside because the knowledge of 



and was found to be some 10^"^ molecules cm ^, see 
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the actual electron energy distribution function and its spatial distribution may be of im- 
portance for a correct extraction of the information about plasma parameters. Langmuir 
probe measurements of electron energy distributions are not a problem nowadays, even in 
the cases of space-charge layers [Q, RF and microwave discharges (see p4| and 



Refs. therein). We limited ourselves with intensity measurements in the more or less easy 
detectable near-UV region because the following second restriction of our attempt seems to 
be more significant. 

2) Our analysis of the known characteristics of radiative and collisional processes involved 
into the mechanism of formation of observable continuum emission led us to the conclusion 
that in spite of numerous efforts the situation is still unsatisfactory. Therefore, only simple 
models can be used. The corona-like model of the excitation, neglecting the rotational 
structure of upper vibronic levels, and the Franck-Condon approximation for electron impact 
excitation are actually available and justified right now. 

Nevertheless, even in such complicated conditions we were happy to find that more or 
less accurate measurements of relative and absolute spectral distributions of the continuum 
intensity may be used for certain estimations of the vibrational temperature, for the ana- 
lysis of the excitation mechanisms and for getting information about the rate of hydrogen 
dissociation by electron impact. 

The last one looks more prospective because in some cases it can be already used as it is. 
If, for example, one is interested in the relative dependence of the electron impact dissociation 
rate on the position in plasma or on the discharge conditions, the information may be easily 
obtained just by a measurement of the continuum emission intensity even without any 
calibrations. The method of the determination of the absolute rate of the electron impact 
dissociation may be considerably improved, if the electron energy distribution function would 
be determined experimentally or calculated somehow. Then there would be no need to use 
our graph in Fig. but to calculate the values in an explicit way by Eq. @. 

This rather easy and independent method of determination of the absolute values of the 
rate of electron impact dissociation of hydrogen may be used in gas discharge physics and 
in plasma chemistry within two approaches. 

1) It may be used for the determination of the hydrogen dissociation degree (which is 
often considered as a very important plasma parameter), if one is able to determine the 
effective mean lifetime of atoms in plasma with certain theoretical model of diffusion and 
association of atoms (see e.g. [|65[] ). 

2) If the value of the dissociation degree is determined by one of the existing methods 
(see |6^j6^) then the effective lifetime of atoms in plasma may be determined with the help 
of the rate of dissociation and the corresponding balance equation for atoms.Q 

In both cases the method proposed in the present work is giving new opportunities. 
Further improvements of the methods developed in present work need more detailed 



Laser-induced fluorescence spectroscopy using two-photon excitation is now being planned to 
determine the concentration of hydrogen atoms under the same conditions in plasma. In combina- 
tion with the results derived in present paper by emission spectroscopy new channel of information 
about the lifetimes of hydrogen atoms in plasma can be opened. 
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information about the elementary processes especially those about the cross sections of 
electronic-vibrational (better electronic-vibro-rotational) a^E+,f', A^' <— X^H'^ ,v, N elec- 
tron impact excitation and the rate coefficients of the quenching of a^T^^ ,v' , N' levels in 
collisions with electrons, atoms and molecules. We certainly understand those are not sim- 
ple problems, but we may hope that the needs of plasma diagnostics and our present work 
may be considered as some kind of the challenge for our colleagues working in the field of 
electronic and atomic collision physics. 
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FIG. 1. Spectral distribution of a^S+, —i- 6^5]+ spontaneous emission transition probabilities 
for various initial vibronic levels f ' = — 3 (curves — 3) of the H2 molecule ||6|,^ . 
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FIG. 2. Relative distributions of the continuum intensity measured in the capillary- arc dis- 
charge with / = 100 mA (curve 1) and those calculated for Tvib=0, 3000, 5000 K (curves 2, 3, and 
4). 
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FIG. 3. Normalized continuum intensity J^^(A, Ao=355.5 nm, T^ib) calculated by Eq. (|l5|) for 
various vibrational temperatures Tvib = 0, 2000, 3000, 4000, 6000 K (curves 1, 2, 3, 4, and 5) and 
measured in capillary arc discharge at currents i = 50, 200, and 300 mA (curves 6, 7, and 8). 
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FIG. 4. Normalized continuum intensity J^^(A, Ao=355.5 nm, Tvib) calculated for A = 300 nm 
as a function of the vibrational temperature Tvib. 
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FIG. 5. Experimental records (counts of CCD detector) of the continuum emission obtained 
in the dc capillary-arc discharges in (a) D2, (b) H2 with variation of the discharge current i = 50, 
100, 200, 300, 400, and 500 mA (curves 1, 2, 3, 4, 5, and 6) and (c) in the microwave discharge 
{p = 0.5 mbar, H2:Ar=l:l). Spectral distribution of the sensitivity is not taken into account. The 
measurements have been performed with gratings 600 (I) and 1200 (II) grooves per mm. 
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FIG. 6. Signals of the continuum emission for Aq = 300 and 355.5 nm (points 1 and 2) against 
discharge current of H2 capillary arc. 




i (mA) 

FIG. 7. Vibrational (1,2) and rotational (3,4) temperatures in the ground X^'^'^ state of the 
H2 in the capillary-arc pure hydrogen plasma for various discharge currents i determined from the 
continuum intensity distribution (1) and from the emission of Fulcher-a bands (2-4). The data 
(1-3) - present work, (4) - Ref. ||. 
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FIG. 8. The relative continuum intensity (normalized for unity at A = 400 nm) measured in a 
H2+Ar microwave plasma for two different ratios of the components at ^? = 0.5 mbar. 




p (mbar) 

FIG. 9. The scaled population densities of ^'^P} levels of Ar as a function of total pressure 
p for constant gas mixture H2+Ar (1:1) in the microwave discharge plasma. The total average 
population of all 3s^3p^4s levels may be obtained by multiplication of 12 to the dashed curve. 
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FIG. 10. Spectral distribution of the continuum emission in the H2+Ar (1:1) microwave plasma 
for the total pressure p = 0.5 mbar. 1 - experimental data points, 2 - calculated by Eq. (|2^ ) with 
the optimal set of parameters A, ci,C2,C3 shown in Table IX, 3 - the contribution of the electron 
impact excitation for Tvib — > with v'^^-^='i, 4 - contribution of the Ar*^H2 excitation transfer 
to a?T,+,v'=Q. 
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FIG. 11. Ratios of the rate coefficients of electron impact excitation of b^T^u {(^b^x)-, ^^^g^ 
{aa^x) states and of the total electron impact dissociation (odiss) calculated as a function of the 
electron temperature. 
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FIG. 12. Radiative dissociation rate (d[ii2]/dt)ab (D) and total electron impact dissociation 
rate ((i[H2]/(ii)totai (•) obtained in the H2+Ar microwave discharge: (a) as a function of the total 
pressure p for constant gas mixture (1:1) and (b) for various gas mixtures under constant total 
pressure p = 0.5 mbar. 
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TABLES 



TABLE L Spectral distributions of transition probability for =0-6^ b'^Ti'^ sponta- 

neous emission of the H2 molecule and relative continuum intensity calculated for Tyib — > by 
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TABLE II. 
probabilities of d'^Il~,v, N- 



Radiative lifetimes of d^H^ , v, N= 



-1 levels and spontaneous emission transition 



=1 — > a^Ti'^ ,v, N=l spectral lines of H2. 



Tdvi (ns) 



present work'^ 



present work 




1 
2 
3 
4 
5 
6 



40.7 ± 1.4 

38.4 ± 1.3 

39.5 ± 1.9 
39.5 ±0.9 
19.0 ± 1.0 
15.2 ± 1.2 
16.0 ±3.0 



38.7 ±2.0 
39.7 ±2.1 
40.9 ±2.2 
42.2 ±2.2 
43.6 ± 2.4 
45.5 ±2.5 



24.4 ± 1.4 
20.6 ±1.1 
16.9 ± 1.0 
13.6 ±0.8 
10.6 ±0.6 
8.1 ±0.5 
5.9 ±0.4 



'^Experiment. 

^Semiempirical calculation in the adiabatic approximation with dipole moment from |1C]. 



TABLE III. Experimental and ab initio data on the radiative lifetimes of levels of 

H2, p. w. - present work. 
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12.4 
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10.3 



9.29 



8.64 



..07 



7.63 



p.w. 



experimental 



11.0 ±0.42 

26 ±2 
11.9 ± 1.2 



9.94 ±0.39 



10.6 ±0.6 
10.8 ± 1.1 



9.1 ± 1.0 



35±8 



10±2 

. 10.45±0.25 



9.62±0.20 



TABLE IV. Relative cross sections for electron impact excitation (in the maximum) of the 
A^'=l levels from X^Tj^ ,v=0,N=l level and corresponding ratios of Franck-Condon fac- 



d'U-,v' 
tors. 



v' 






'^xoi/^xdi 






FCF 


|39| 


|4C] 


132] 


i 


present work 


1 


42| 





0.57 


0.68 


0.52±0.07 


0.66±0.07 


0.70±0.07 


0.55 


1 


0.86 


1.09 


0.91±0.13 


0.88±0.09 


0.96±0.09 


0.95 


2 


1.00 


1.00 


1.00±0.14 


1.00±0.10 


1.00±0.10 


1.00 


3 


0.78 


0.86 


0.96±0.14 


0.77±0.08 


0.84±0.08 


0.86 


4 








0.35±0.04 


0.87±0.09 


0.66 


5 








0.19±0.02 


0.60±0.06 


0.47 


6 








0.11±0.01 


0.32±0.03 


0.33 
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TABLE V. Franck-Condon factors for a^Ej, v' ^ X^T.+,v transitions of Ha [||]. 



















v' 


= 


1 


2 


3 


4 


5 


■u = 6 





0.20761 


0.39958 


0.28411 


0.09341 


0.01439 


0.00089 


0.00001 


1 


0.25478 


0.06483 


0.08307 


0.32592 


0.21741 


0.05018 


0.00378 


2 


0.20249 


0.00503 


0.16438 


0.00281 


0.20634 


0.30686 


0.10324 


3 


0.13470 


0.06063 


0.05177 


0.08235 


0.07231 


0.07873 


0.34105 


4 


0.08235 


0.09580 


0.00050 


0.10148 


0.00727 


0.12529 


0.01180 


5 


0.04840 


0.09603 


0.01522 


0.04193 


0.06678 


0.00933 


0.11680 


6 


0.02803 


0.07908 


0.04137 


0.00439 


0.07005 


0.01583 


0.04615 



TABLE VL Franck-Condon factors for d^B^, t>' ^ i; transitions of H2 |2). 



v' 
















v = a 


1 


2 


3 


4 


5 


V = 6 





0.09995 


0.29272 


0.33797 


0.19710 


0.06160 


0.00994 


0.00070 


1 


0.17029 


0.16255 


0.00036 


0.16128 


0.29619 


0.16752 


0.03860 


2 


0.18029 


0.02851 


0.08310 


0.10636 


0.00914 


0.23640 


0.26345 


3 


0.15451 


0.00075 


0.10888 


0.00041 


0.12133 


0.02528 


0.11911 


4 


0.11833 


0.02432 


0.05755 


0.03636 


0.05553 


0.04245 


0.09099 


5 


0.08505 


0.05077 


0.01334 


0.06851 


0.00096 


0.08332 


0.00043 


6 


0.05899 


0.06393 


0.00002 


0.05697 


0.01538 


0.03661 


0.04690 



TABLE VIL Results of x^-minimisation of the experimental data shown on Fig. ^ with various 
number of adjusted parameters (t'max+1)- 



X P 



( 



max 



57.1 10-^ 0.093 

1 9.89 0.62 0.108 

2 2.10 0.73 0.124 

3 0.426 0.84 0.127 

4 0.308 0.98 0.833 

5 0.149 1.00 6.670 

6 0.159 1.00 2.820 
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TABLE VIII. Experimental and calculated with v'^^^ = 3 relative population densities for 
^vib — > of the a^T,'^,v' vibronic states and relative intensities [Eq. (|l^)] in the range 
A=230-400 nm. Pure hydrogen capillary-arc discharge i=100 mA - see Fig. ||. Error(s) in brackets 



correspond to last digit (s). 



V 




AT 1 AT 














expt. 




calc. 




expt. 




calc. 





L00(3) 




1.00 




0.50(2) 




0.50 


1 


L20(4) 




1.14 




0.27(1) 




0.25 


2 


0.79(7) 




0.81 




0.14(1) 




0.14 


3 


0.86(11) 




0.52 




0.09(1) 




0.06 


4 






0.28 








0.03 


5 






0.15 








0.01 


6 






0.08 








0.01 


TABLE IX. 


Relative contributions 


Ci of various 


channels of the excitation of the continuum 


emission [Eq. 


I)] in the microwave plasma with constant pressure p = 0.5 mbar 


for various H2+Ar 


mixtures and the rates of radiative dissociation —{d\H2]/dt)ab 


in 10^^ cm ^ s ^ 


Error 


in brackets 


corresponds to last digit. 














Model 


[Ar]:[H2] 


Cl 


C2 






-(d[H2]/dt)a6 














(1)- 


(2)'^ 




1:4 


0.90(4) 


0.10(6) 


0.00(3) 


0.93 


1.40 


3.4 




1:1 


0.79(4) 


0.21(6) 


0.00(3) 


1.04 


1.05 


2.4 


ci, C2, C3 7^ 


4:1 


0.74(5) 


0.23(7) 


0.03(3) 


1.27 


0.81 


1.8 




9:1 


0.72(5) 


0.25(7) 


0.03(4) 


1.66 


0.78 


1.7 




19:1 


0.77(6) 


0.18(8) 


0.05(4) 


1.90 


0.70 


1.6 




1:4 


0.90(4) 


0.10(4) 




0.89 


1.37 


3.3 


ci, C2 ^0 


1:1 


0.79(4) 


0.21(4) 




0.98 


1.05 


2.4 


and 


4:1 


0.73(5) 


0.27(5) 




1.26 


0.82 


1.8 


C3 = 


9:1 


0.71(5) 


0.29(5) 




1.64 


0.78 


1.7 




19:1 


0.75(6) 


0.25(6) 




1.96 


0.70 


1.6 




1:4 


0.97(3) 




0.03(2) 


1.09 


1.36 


3.4 


ci , C3 ^0 


1:1 


0.93(3) 




0.07(3) 


1.68 


1.06 


2.6 


and 


4:1 


0.90(4) 




0.10(3) 


2.10 


0.82 


2.0 


C2 = 


9:1 


0.89(4) 




0.11(3) 


2.58 


0.79 


1.9 




19:1 


0.90(4) 




0.10(3) 


2.35 


0.71 


1.7 



'^Integrated over A = 225 — 400 nm. 



Total radiative dissociation rate with extrapolation for entire range of A. 
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